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Improved search for muon-neutrino to electron-neutrino oscillations in MINOS 
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We report the results of a search for appearance in a Vfi beam in the MINOS long-baseline 
neutrino experiment. With an improved analysis and an increased exposure of 8.2 x 10^° pro- 
tons on the NuMI target at Fermilab, we find that 2 sin2(6l23) sin^(26li3) < 0.12 (0.20) at 90% 
confidence level for 5=0 and the normal (inverted) neutrino mass hierarchy, with a best fit of 
2sin2(6'23)sin^(26li3) = 0.04lj:[];J!3i (O.OTQto.oss)- The 013=0 hypothesis is disfavored by the MINOS 
data at the 89% confidence level. 



PACS numbers: 14.60.Pq, 14.60.Lm, 29.27.-a 



2 



It has been experimentally established that neutrinos 
undergo flavor change as they propagate [IHZj • This phe- 
nomenon is well-described by three-flavor neutrino oscil- 
lations, characterized by the spectrum of neutrino masses 
together with the elements of the PMNS mixing ma- 
trix [5| . This matrix is often parametrized by three Euler 
angles 9ij and a CP-violating phase 5. While 9i2 and ^23 
are known to be large PQIHH], ^13 appears to be relatively 
small [9Ul3|. with the tightest limits so far coming from 
the CHOOZ Ho] and MINOS [12] experiments. The T2K 
collaboration has recently reported indications of a non- 
zero value for ^13 at the 2.5(7 confidence level (C.L.) [H] . 
This letter reports new ^13 constraints from the MINOS 
experiment, using an increased data set and significant 
improvements to the analysis. 

MINOS is a two-detector long-baseline neutrino os- 
cillation experiment situated along the NuMI neutrino 
beamline il5\ The 0.98-kton Near Detector (ND) is lo- 
cated on-site at Fermilab, 1.04 km downstream of the 
NuMI target. The 5.4-kton Far Detector (FD) is located 
735 km downstream in the Soudan Underground Labo- 
ratory. The two detectors have nearly identical designs, 
each consisting of alternating layers of steel (2.54 cm 
thick) and plastic scintillator (1 cm). The scintillator lay- 
ers are constructed from optically isolated, 4.1 cm wide 
strips that serve as the active elements of the detectors. 
The strips are read out via optical fibers and multi-anode 
photomultiplier tubes. Details can be found in Ref. [I6] . 

The data used in this analysis come from an expo- 
sure of 8.2x10^" protons on the NuMI target. The cor- 
responding neutrino events in the ND have an energy 
spectrum that peaks at 3 GeV and a flavor composition 
of 91.7% v^, 7.0% V^, and 1.3% v^+^e, as estimated 
by beamline and detector Monte Carlo (MC) simula- 
tions, with additional constraints from MINOS ND data 
and external measurements [T7] . The two-detector ar- 
rangement and the relatively small intrinsic compo- 
nent make this analysis rather insensitive to beam uncer- 
tainties. Neutrino-nucleus and final-state interactions are 
simulated using NEUGEN3 [TH], and particle propagation 
and detector response are simulated with GEANT3 [19 . 



MINOS is sensitive to 6'i3 through Vu 



Vp oscilla- 



tions. To leading order, the probability for this oscilla- 
tion mode is given by 



sin2(6l23) sin2(2e'i3) %va^(\.TJ l\m\^LlE) , 



where Ato^j (in units of eV^) and ^23 are the domi- 
nant atmospheric oscillation parameters, L (in km) is the 
distance between the neutrino production and detection 
points, and E (in GeV) is the neutrino energy. We set 
constraints on ^13 by searching for an excess of Vf^ events 
at the FD. Matter effects and possible leptonic CP vi- 
olation modify the above probability significantly [20] . 
hence our results are presented as a function of b and the 
neutrino mass hierarchy. 



Events in the MINOS detectors can be characterized 
by the spatial patterns of energy deposition in the scin- 
tillator strips. Charged-current (CC) interactions are 
identified by a muon track extending beyond the more lo- 
calized hadronic recoil system. A single detector plane is 
1.4 radiation lengths thick, so the electron from a CC 
interaction penetrates only a few planes (typically 6—12), 
leaving a transversely compact pattern of activity inter- 
mingled with the associated hadronic shower. Neutral- 
current (NC) interactions can mimic this pattern, partic- 
ularly when neutral pions are present. 

To obtain a z^e-enriched sample, we apply a series of 
selection criteria to the recorded neutrino events. We 
require that the neutrino interaction occur within a fidu- 
cial volume. We eliminate most CG interactions by 
rejecting events with a track longer than 24 planes and 
events with a track extending more than 15 planes be- 
yond the hadronic shower. We require that an event have 
at least five contiguous planes with energy greater than 
half that deposited on average by a minimum ionizing 
particle. The calorimetrically determined event energy 
must lie between 1 and 8 GeV, as events below 1 GeV are 
overwhelmingly from NC interactions and events above 
8 GeV have negligible — )■ Vf. oscillation probability. 
The time and reconstructed direction of each event must 
be consistent with the low-duty-cycle NuMI neutrino 
source. These "pre-selection" criteria preserve 77% of 
oscillation-induced CC events originating in the fidu- 
cial volume while passing 8.5% of z/^ CC, 39% of NC, 
54% of v-r CC, and 35% of intrinsic v,, CC events, as 
estimated by the simulation. 

Further background suppression requires a more so- 
phisticated examination of the energy deposition pat- 
terns. Earlier MINOS appearance searches used an 
artificial neural network event classifier with eleven input 
variables characterizing the transverse and longitudinal 
profiles of an event's activity in the detector [TT ] fT2 l [2T ] . 
The present analysis uses a nearest-neighbors algorithm 
dubbed "library event matching" (LEM) [Hj. In LEM, 
each candidate event is compared to 5 x 10^ simulated sig- 
nal and background events [53] , one by one, to find the 50 
that look most similar to the candidate event. A library 
event is rejected if its reconstructed energy, number of 
active strips, or number of active planes differs from that 
of the candidate event by more than 20%. The similarity 
of the candidate event to each remaining library event is 
quantified by the following likelihood: 



log£ 



Natrips 

E 



log 



P«,„d;A)P«b;A)dA 



where the sum is taken over all strips with a signal above 
3 photoelectrons in either of the two events, n\ is the 
charge (in photoelectrons) observed on strip i in event 
X (with X either "candidate" or "library"), and P{n;\) 
is the Poisson probability for observing n given mean A. 
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Since events occur throughout the detector volume, each 
event is translated to a fixed reference location before 
£ is evaluated. Strips far away from the event's central 
axis are combined before comparison. Additionally, li- 
brary events are shifted by ±1 plane in search of a better 
likelihood. 

The final classifier is formed using a neural network 
that takes as its inputs the reconstructed event energy 
along with three variables derived from the best-match 
ensemble: (1) the fraction of the 50 best-matched events 
that are true v,. CC events, (2) the average inelasticity 
{y) of those v^, CC events, and (3) the average fraction 
of charge that overlaps between the input event and each 
j'e CC event. The resulting LEM discriminant is shown 
in Fig. [1} 



Uncertainty source 


Uncertainty on 
background events 


Event energy scale 


4.0% 


u-r background 


2.1% 


Relative FD/ND rate 


1.9% 


Hadronic shower model 


1.1% 


All others 


2.0% 


Total 


5.4% 



TABLE I: Systematic uncertainties on the number of pre- 
dicted background events in the FD in the signal region, de- 
fined by LEM>0.7. The final ^13 measurement uses multiple 
LEM and reconstructed energy bins and thus uses a full sys- 
tematics covariance matrix. These uncertainties, which are 
small compared to the statistical errors, lead to a 7.0% loss in 
sensitivity to sin^(2&i3). The "All others" category includes 
uncertainties relating to the neutrino flux, cross sections, de- 
tector modeling, and background decomposition. 
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FIG. 1: [Top] Distribution of the LEM discriminant for events 
in the Near Detector that pass the pre-selection requirements. 
Data (points) and Monte Carlo simulation (histogram) are 
shown, with the magnitude of the systematic uncertainty in- 
dicated by the band. This uncertainty is highly correlated 
between the ND and FD and thus cancels out to a large de- 
gree when we form our FD predictions. [Bottom] Expected 
background and signal distributions in the Far Detector for 
sin^(26i3)=0.1. The signal distribution has been multiplied 
by 10 for visibility. 

We form a prediction for the FD event rate, in each of 



15 bins (specified below) of LEM discriminant and recon- 
structed energy, using the corresponding rate observed in 
the ND. The ND rates are first broken down into indi- 
vidual background contributions, as different background 
types translate differently from the ND to the FD due to 
oscillations and beamline geometry. To determine the rel- 
ative background contributions in the ND rates, we apply 
the Ve selection to ND data collected in multiple beam 
configurations with differing neutrino energy spectra and 
thus differing background compositions. This allows the 
construction of a system of linear equations that can be 
solved for the relative contributions of NC, CC, and 
intrinsic CC backgrounds in the primary low-energy 
beam configuration (12j . The measured composition of 
ND events, averaged over the range LEM>0.7 for recon- 
structed energy between 1 and 8 GeV, is (61 ± 1)% NC, 
(24 ± 1)% CC, and (15 ± 1)% CC. 

We convert the resulting decomposed ND rates directly 
into predictions for the FD rates using a Monte Carlo 
simulation. More specifically, we use the simulated ra- 
tio of FD and ND rates, for each background type and 
for each LEM and energy bin, as the conversion factor 
for translating the measured ND rate into the FD pre- 
diction. We evaluate uncertainties on these ratios using 
systematically modified samples of simulated ND and FD 
data. The dominant systematic effects are summarized 
in Table n 

Since the ND collects negligibly few events arising from 
— or — > Vt oscillations, the FD rates for 
these events are estimated using the simulation plus the 
observed t'^ CC rates in the ND. For v^^ CC events, 
we further apply an energy- and LEM-dependent cor- 
rection to the FD predictions that is derived from hy- 
brid events composed of electrons from simulation and 
hadronic showers from data. The hadronic showers are 
obtained by removing the muon hits from cleanly iden- 
tified CC events [HI US] , and the electromagnetic 
shower simulation is verified using a pure sample of elec- 



Event class 



sm"(26li3) 
0.1 



MINOS Far Detector Data 



NC 


34.1 


34.1 


CC 


6.7 


6.7 


V, CC 


6.4 


6.2 


Vr CC 


2.2 


2.1 


VlL -> J'e CC 


0.2 


19.1 


Total 


49.6 


68.2 



TABLE II: Expected FD event counts for LEM>0.7, assum- 
ing 6*23 =f, Amia = 2.32x10"^ eV^ and 5 = 0. The first 
line refers to the intrinsic component in the beam. In the 
^13=0 case, a small amount of — !> oscillation occurs due 
to non-zero Am^i. 



trons recorded by the MINOS Calibration Detector [55]. 
The breakdown of expected FD events is given in Ta- 
ble [n] An analysis of beam-off detector activity yielded 
no z^e candidate events, resulting in a 90% C.L. upper 
limit on cosmogenic backgrounds in the primary analysis 
region of 0.3 events. We find that (40.4 ± 2.8)% of CC 
signal events end up in the signal region, LEM>0.7. 

Most of the analysis procedures can be tested directly 
on two signal-free or near-signal-free sideband samples. 
First, the "muon-removed" hadronic showers described 
above, before they are merged with simulated electrons, 
represent a sample of NC-like events. The predicted 
and observed LEM distributions in the FD agree for 
this sample, with x^/-^d.o.f. =9.7/8 using statistical errors 
only. Second, FD events satisfying 0<LEM<0.5 make 
up a background-dominated sample for which we predict 
370 ± 19 background events (statistical error only). We 
observe 377 events, in agreement with prediction. Form- 
ing the prediction for the latter sideband exercises all 
aspects of the analysis up to the final signal extraction, 
including the full ND decomposition procedure and the 
ND-to-FD ratios derived from simulation. 

In previous MINOS analyses [UlIE], the v^^ appear- 
ance search was conducted by comparing the total num- 
ber of J'e candidate events in the FD to the expected 
background. A similar approach applied to the present 
data yields 62 events in the signal region of LEM>0.7, 
with an expectation of 49.6 ± 7.0(stat.) ± 2.7(syst.) if 
013=0. However, we gain 12% in sensitivity by fit- 
ting the FD sample's LEM and reconstructed energy 
(^-roco) distribution in 3x5 bins spanning LEM>0.6 and 
1 GeV < i?rcco < 8 GeV. The energy resolutions for 
hadronic and electromagnetic showers at 3 GeV are 32% 
and 12%, respectively 16]. Figure [2] shows the FD data 
and predictions used in the fit, along with the extracted 
best-fit signal. 

Figure [3] shows the regions of oscillation parame- 
ter space allowed by these data. For the fit, we 
use a three-flavor oscillation framework [5D] includ- 
ing matter effects [17j, and we use the Feldman- 
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FIG. 2: Reconstructed energy spectra for i^e CC candidate 
events in the Far Detector. The black points indicate the 
data with statistical error bars shown. The histogram indi- 
cates the expected background (unfilled area) together with 
the contribution of i^^ — > signal (hatched area) for the 
best-fit value of sin2(26li3) = 0.041. 



Cousins procedure 
gions. We assume 



to 



re- 



|Am§2| 



calculate the allowed 
= (2.32;0;J^)xl0-3 eV^ [6], 
(7-59lo;2?)xl0-5 eV^ [T], 6(23 = 0.785 ±0.100 g], 
and 6*12 = 0.60 ±0.02 [I]. The influence of these oscilla- 
tion parameter uncertainties is included when construct- 
ing the contours. 

Prior to unblinding the FD data, we planned to fit only 
the LEM distribution integrated over energy. However, 
the excess over background in the upper energy range 
prompted the inclusion of energy information so that the 
fit could weigh events appropriately when extracting ^13 
constraints. If we had performed the signal extraction 
over LEM bins only, the best fit and 90% C.L. upper 
limit for sin^ (26*13) would each change by +0.006. A thor- 
ough study of high-energy events in the signal and side- 



5 




Am^ > 

— MINOS Best Fit 

68% C.L. 
■ 90%C.L. 
■■■■ CHOOZ90%C.L. 

2s\n\^ = 1 for CHOOZ 



Am^< 



MINOS 
8.2x10^° POT 



0.3 



0.4 



2sin^(2ei3)sin\3 



FIG. 3: Allowed ranges and best fits for 2 sin (e23)sin2(2ei3) 
as a function of S. The upper (lower) panel assumes the nor- 
mal (inverted) neutrino mass hierarchy. The vertical dashed 
line indicates the CHOOZ 90% C.L. upper limit assuming 
023^"^ and Amia = 2.32x10"^ eV^ HU]. 



band samples, including events between 8 and 12 GeV, 
indicates that the high-energy predictions are robust and 
that the selected events are free of irregularities. 

In conclusion, using a fit to Ve discriminant and 
reconstructed energy 2D distribution of FD i/g can- 
didate events, we find that 2 sin^(^23) sin^(26'i3) — 
0.04lt°;°*J (0-079l[!;|^^3) for the normal (inverted) 
mass hierarchy and <5=0. We further find that 
2sin2(6l23)sin2(20i3)<O.12 (0.20) at 90% C.L. Using the 
less sensitive techniques of the 2010 analysis [T^] on the 
current data set yields a consistent measurement |29j . 
The 013=0 hypothesis is disfavored by the MINOS data 
at the 89% C.L. This result significantly constrains the 
range allowed by the T2K data [TJ and is the most 
sensitive measurement of ^13 to date. 
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